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多 策略 增强 花 授 粉 算法 及 其 应 用 
李 大 海 ， 伍 兆 前 ， 王 振东 
(江西 理工 大 学 信息 工程 学 院 , 江西 闽 州 341099) 


摘 要 : 针对 经 典 花 授粉 算法 容易 陷入 局 部 最 优 解 和 收 你 速度 慢 的 缺点 ， 提 出 一 种 增强 型 透镜 成 像 策 略 和 随机 邻 域 
变异 策略 的 花 授 粉 算 法 。 通 过 增强 型 透镜 成 像 策略 扩展 花 授 粉 算 法 的 搜索 空间 ， 增 加 解 的 多 样 性 有 助 于 算法 跳出 局 
部 最 优 解 。 引 入 随机 邻 域 变异 策略 ， 和 借助 邻 域 内 的 信息 指导 算法 搜索 ， 增 强 算法 的 收 人 敏 精度 和 搜索 速度 。 并 且 对 改 
进 后 的 花 授粉 算法 和 4 种 其 他 改进 算法 在 CEC2013 测试 函数 上 比较 ， 实 验证 明 改 进 后 的 多 策略 花 授 粉 算 法 不 管 是 
在 收敛 精度 还 是 搜索 速度 都 比 对 比 算 法 优秀 。 最 后 把 多 策略 花 授粉 算法 应 用 在 汽车 传动 参数 模型 上 研究 该 算法 的 实 
际 效用 ， 结 果 表 明 多 策略 花 授 粉 算 法 在 汽车 传动 参数 优化 问题 上 都 优 于 对 比 工法 。 
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Multi-strategy flower pollination optimization algorithm for vehicle power transmission parameters 


Li Dahai, Wu Zhaogian, Wang Zhendong 
一 一 (School of Information Engineering, Jianexi University of Science & Technology, Ganzhou Jianexi 341099, China) 


Abstract: Classic flower pollination algorithm (FPA) can be easily exposed to the shortcomings of local optimal solution and 
slow convergence velocity. In view of these shortcomings, this paper proposed an FPA with an enhanced lens imaging strategy 
and random neighborhood-based mutation strategy. The lens imaging strategy can help the algorithm avoid the shortcoming 
of local optimal solution by expanding the search space of FPA to increase the diversity of the solution. The introduction of 
random neighborhood-based mutation strategy can enhance the convergence accuracy and search speed of the algorithm by 
guiding algorithm search with information in the neighborhood. A comparison of the improved FPA with four other improved 
algorithms on CEC2013 test function found that the improved multi-strategy FPA performs better than the comparison 
algorithms in both convergence accuracy and search speed. To study its practical utility, this paper applies the multi-strategy 
FPA into the automobile transmission parameter model and the results indicate that multi-strategy FPA is better than the 
comparison algorithm in optimization of automobile transmission parameters. 

Key words: random neighborhood variation; lens imaging; flower pollination algorithm; parameter optimization; 
convergence accuracy 


5 0 引言 自 适 应 的 贪 禁 策 略 ， 其 中 全 局 性 精英 对 立 的 学 习 增 强 了 和 群体 
C= 加 的 多 样 性 ， 局 部 自 适 应 的 贪 焚 策 略 增强 了 其 探索 能 力 。 文 献 
CS 花 授 粉 算法 (Flower Pollination Algorithm, FPA) 中 是 [10] 中 提出 了 一 种 结合 差分 进化 变异 策略 的 ssFPA/DE 算法 ， 
Yang 于 2012 年 提出 一 种 简单 高 效 的 元 启发 式 群 智 能 优化 算 ”该 新 算法 将 FPA 算法 的 搜索 策略 特征 和 差分 进化 的 变异 策 
法 。 该 算法 将 植物 的 异 花 授 粉 和 自 花 授粉 的 过 程 类 比 到 算法 略 相 结合 从 而 增强 算法 的 搜索 效率 。 文 献 [11] 中 针对 FPA 算 
中 的 全 局 搜索 和 局 部 搜索 ， 平 衡 全 局 搜索 和 局 部 搜索 之 间 的 。 ”法 易 陷 入 局 部 最 优 ， 后 期 收敛 速度 慢 的 缺点 ， 提 出 一 种 基于 
切换 。 由 于 FPA 算法 具有 结构 简单 、 参 数 少 、 鲁 棒 性 和 适应 引力 搜索 机 制 的 GSFPA 算法 ， 通 过 花 条 个 体 间 的 引力 和 本 
性 强 等 特点 ， 其 已 经 应 用 于 求解 诸多 领域 的 复杂 优化 问题 求 。 ” 身 算 法 的 Levy 飞 行 共同 实现 个 体位 置 的 更 新 , 使 花灯 受 Levy 
解 ， 如 无 线 传 感 飞行 和 个 体 间 引力 的 双重 影响 ,个体 之 间 通 过 共享 优化 信息 
器 网 络 布置 优化 所 ,混凝土 3D 框架 优化 中 ， 医 学 图 像 分 。 提高 算法 的 寻 优 能 力 。 文献 [12] 中 提出 FA-FPA 算法 , 在 FPA 
制 外 ， 切 比 雪 夫 多 项 式 神经 网 络 优化 外 。 但 FPA 算法 也 有 着 。 算法 的 局 部 授粉 过 程 中 引入 自 适 应 的 变异 因子 , 并 对 FPA 算 
和 其 他 群 智能 优化 算法 相同 的 问题 ， 容 易 陷入 局 部 最 优 解 并 ”法 的 概率 转换 进行 自 适应 调整 。 最 后 将 其 与 萤火虫 算法 相 结 
收敛 速度 慢 四 ， 国 内 外 学 者 针对 FPA 算法 已 经 提出 诸多 的 。” 合 ， 实 验 表明 FA-FPA 算法 具有 更 高 的 收敛 性 和 稳定 性 。 
改进 措施 ,文献 [7] 提 出 一 种 结合 克隆 选择 策略 的 改进 FPA 入 上 述 对 于 FPA 改进 的 算法 主要 外 策 
> 


dk 


主要 集中 在 三 点 : 一 是 借用 策 
法 其 使 用 随机 游 走 蔡 代 FPA 算法 中 使 用 的 Levy 飞行 并 且 在 略 某 方面 的 优越 性 弥补 FPA 算法 的 自身 缺陷 。 例如 模拟 退火 
局 部 授粉 之 前 克隆 最 优 解 ， 以 提高 最 优 解 的 接受 概率 。 文 献 的 突 跳 机 制 使 其 避免 陷入 局 部 最 优 解 ， 引 力 搜索 机 制 调整 个 
[8] 提出 了 一 种 结合 模拟 退火 策略 的 SFPA 算法 , 该 增强 FPA 体位 置 更 新 。 二 是 增加 种 群 多 样 性 。 例 如 全 局 性 精英 对 立 学 
算法 利用 模拟 退火 策略 使 算法 避免 陷入 局 部 最 优 解 ， 从 而 提 习 增 强 种 群 多 样 性 。 三 是 融合 其 他 算法 ， 借 用 其 他 算法 的 优 
高 了 算法 的 全 局 搜索 能 力 。 文 献 [9] 中 提出 了 一 种 基于 精英 对 秀 机 制 增强 FPA 算法 。 
立 策略 的 EOFPA 算法 该 算法 采用 了 全 局 性 精英 对 立 和 局 部 尽管 目前 已 经 提出 了 诸多 的 改进 FPA 算法 ， 但 FPA 算 
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法 仍 存在 诸多 的 改进 空间 。 本 文 提出 一 种 基于 多 策略 的 


FPA 算法 ， 其 主要 从 以 下 两 方 


而 改进 FPA 算法 : 


a) 通过 增强 型 透镜 成 像 策 


各 去 增强 算法 的 种 群 多 样 性 。 


透镜 成 像 策 略 是 一 种 基于 光学 原理 的 反 
能 够 增强 种 群 的 多 样 性 改善 算法 的 搜索 


局 部 最 优 解 。 本 文 提出 一 种 增强 型 透镜 成 像 策 
民 为 一 回合 。 在 标准 透镜 成 像 生成 反 向 


向 学 习 策略 ， 该 策略 
目 ， 帮 助 算法 跳出 
各 ， 以 10 次 迭 
解 的 基础 上 ， 再 利用 


沁 团 ， 


算法 生成 的 最 优 解 和 透镜 成 像 的 反 向 解 


生成 精英 集合 ， 对 精英 集合 内 部 的 


序 对 分 别 求 适应 值 ， 精 


两 者 之 间 的 笛 卡尔 积 


英 集合 内 部 始终 保持 10 个 个 体 ， 新 的 优质 解 会 蔡 代 劣质 解 。 


能 够 平衡 全 局 和 局 部 搜索 ， 
略 和 FPA 粉 算 法 的 搜索 策 
法 的 全 局 搜索 和 局 
适应 ， 减 
邻 域 集合 内 部 的 


CEC2013 测试 函数 上 的 性 


b) 引入 差分 进化 算法 区 


所 以 本 文 提 


邻 域 变异 策 


各 。 由 于 邻 域 变 异 策 


各 相 结合 ， 


小 不 合适 的 变异 因子 对 算法 的 
言 息 指导 搜索 ， 加 快 收 
通过 对 比 FFPA，TMFPA，HLFPA 


并 把 算法 应 用 到 汽车 传动 系统 参数 优化 
步 验证 算法 优化 复杂 问题 的 能 


1 


花 授粉 算法 概述 (FPA) 
然 界 花 授粉 的 过 程 


FPA 算法 借鉴 


的 自 花 授粉 和 异 花 授粉 的 行为 。 授 粉 的 媒介 分 为 两 种 ， 生 


者 昆虫 充当 载体 ， 可 在 大 范围 内 传播 ; 


花 
物 和 非 生 物 。 生 物 授 粉 的 过 程 视 为 异 花 授粉 ， 一 般 
者 
为 


部 搜索 的 随机 切换 ， 将 变异 因子 F 改 为 自 


能 ， 验 证 本 文 提出 算法 的 有 效 1 


花 授粉 ， 通常 是 借助 风 实现 传播 。 为 了 更 好 利 


了 适应 FPA 算 


影响 程度 ， 最 后 借用 
敛 速度 和 收敛 精度 。 
，ISSA 四 种 算法 在 
生 ， 
的 工程 问题 上 ， 进 


， 基 本 模拟 了 自然 界 
动物 或 

非 生物 的 自 花 授粉 视 
j FPA 算 
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FPA 算法 作为 一 种 群 智能 优化 算法 ， 同 样 临 如 何平 衡 算 


法 种 群 多 样 性 、 收 敛 速度 和 收敛 精度 的 问题 。 


本 文 提出 利用 


增强 型 透镜 成 像 反 向 学 习 策略 和 基于 随机 邻 域 变异 策略 两 项 


措施 改进 FPA 算法 ， 以 进 
2.1 增强 型 透镜 成 像 反 向 学 习 策 略 


步 加 强 算法 的 收敛 速度 和 精度 。 


凸透镜 成 像 规律 是 一 种 光学 定律 。 如 图 


1 所 示 ， 假 设 物 


体 P, 在 * 轴 上 的 投影 为 x, 高 度 为 h， 通过 凸透镜 可 以 获得 
一 个 实 像 P', 设 实 像 疡 在 x* 轴 上 的 投影 为 x*, 设 投影 * 和 x* 


在 x* 轴 上 的 投影 落 于 区 间 
示 的 关系 : 


(a+b)/2—x _h 
x*—(a+b)/2 h* 


[4 如 之 内 , 则 x 和 x* 满 足 如 式 (4) 所 


(4) 


令 n= 了 好 ， 并 称 n 为 缩放 因子 。 对 式 (4) 进 行 变 换 即 可 得 


h* 
到 x* 的 计算 式 (5): 


水 


_a+D QT+D x 
2 2n n 


若 将 x* 作 为 
x* 可 以 视 为 解 * 的 反 疝 解 。 特 别 的 , 当 n=1 旧 
化 成 式 (6) 的 形式 ， 其 就 是 
X=a+Db—x 
在 透镜 成 像 中 可 以 通过 调整 缩放 因子 nn 站 
某 个 解 的 反 向 解 ， 对 一 般 的 DD 维 优 化 问题 ， 


AS 


(5) 


E 优 化 问题 的 某 个 解 (或 者 解 的 坐标 )， 则 


时 将 ， 式 (5) 可 简 


般 的 反 向 学 习 策略 。 


(6) 
的 大 小 来 动态 获 


设 D 维 向 量 x 


该 优化 问题 的 某 个 解 (或 者 解 的 坐标 )。 


法 解决 优化 问题 , Yang 在 算法 中 假定 每 


生物 充当 载体 ， 并 遵循 Levy 飞行 的 规律 来 实现 ; 
花 授粉 对 应 算法 的 开发 阶段 ， 
同 种 类 的 植物 不 同 花 灯 之 间 通 过 


二 


授粉 ) 为 FPA 算法 的 核心 。FPA 算法 的 全 


世 d 
< 


末 花 只 有 


生成 比例 ; 
4) 概率 转换 参数 pe[0,1] 对 FPA 算法 的 探索 (全 局 授粉 ) 


天 
风 等 因素 影响 ， 整 个 授粉 过 程 中 ， 
上 述 的 描述 可 知 ， 异 花 授 粉 (全 局 授粉 ) 和 自 花 授 粉 (局 部 


此 外 还 需 假定 以 下 4 条 规定 : 


一 株 植物 只 开 一 茶花 ， 


个 花粉 胚 子 ， 一 个 胚 子 对 应 优化 问题 中 的 一 个 


1) 植物 的 异 花 授粉 对 应 算法 的 探索 阶段 , 全 局 搜索 通过 


2) 非 生 物 的 


3) 繁衍 概率 对 应 花 的 特性 ， 两 


F 发 (局 部 授粉 ) 之 间 的 相互 转换 进行 


分 量 D 维 向 量 x* 作为 解 < 


»_Wtb wtb 
认 = 十 一 二 
2 2n n 


FP， % 沁 为 搜索 空间 在 第 j 维 的 上 下 界 。 


I 
ml 


局 部 搜索 是 


风 实 现 传粉 ; 
茶花 (个 体 ) 的 相似 和 关联 


周 节 。 由 于 受到 位 置 、 


局 部 授粉 的 概率 更 大 。 


综合 以 上 的 透镜 成 像 原理 ， 本 文 提 出 一 升 
成 像 策略 。 首 先 将 开始 
5 集合 中 ， 将 其 
并 生成 合 5, 和 5 

Q=S5xS,={(%,%)|h eS and » eS, 


10 次 迭代 生成 的 当前 最 优 解 保存 在 
透镜 成 像 生成 的 反 向 解 保存 在 集合 9 中 。 
的 的 笛 卡 尔 乘积 集合 2 ， 如 式 (8) 所 示 。 


设 态 为 当 向 量 x 的 
分 量 ， 按 式 (7) 得 * 的 反 向 分 量 号 ， 则 可 得 以 冯 作为 的 
的 透镜 反 向 解 。 


(7) 


中 增强 型 的 透镜 


} (8) 


在 随后 的 算法 迭代 中 ， 如 果 当 前 和 帮 代 次 数 不 是 10 的 倍 
数 , 则 将 5，5, 记录 的 个 体 的 适应 值 分 别 和 当前 最 优 解 及 其 反 


J 


X11 =X +yLx(x! —g),rand 


局 授粉 按 式 (1) 进 行 : 
> 万 (1) 


AT(4) sin() 1 
L 一 一 
A Slt4 


其 中 ，x 交 ， 半 为 第 #+1 和 + 代 的 解 ，8 为 全 局 最 优 解 ，7 是 步 
长 控制 因子 , L 是 Levy 飞行 的 主 步 长 5 


步 长 计算 公式 如 下 : 


2) 


>0 


Pp ,44=3/2, TW 是 标准 的 伽 马 函数 。 


于 按 式 (1) 更 新 总 ， 


所 以 整个 全 局 搜索 阶段 都 是 朝 着 当前 和 
使 用 Levy 飞行 每 一 次 都 随机 调整 步 长 , 以 确保 整 
fF 性 。 通 过 此 种 方式 获取 的 每 一 代 的 最 


向 探索 ， 
个 种 群 的 多 刷 
当下 一 次 和 


的 全 局 最 优 解 8 的 方 


CL 解 8 充 


尺 的 基 和 向量， 保证 最 优 解 基 
FPA 的 局 部 授粉 按 如 下 式 (3) 进 行 : 


其 中 ，2 是 在 [0,1] 区 间 上 服从 均匀 分 布 


与 # 不 同 的 随机 选取 的 花粉 。 其 


阶段 日 用 * 扰动 确保 每 次 获取 的 解 都 具 


=X+EX(X—M),rand <p 


因 的 传承 。 


GB3) 
的 随机 数 ， 已 乓 代表 


的 是 


有 较 好 的 开发 性 ， 


有 随机 性 。 


确保 算法 在 局 部 搜索 


再 按 式 (8) 重 新 生成 集合 2 。 


向 解 的 适应 值 进行 对 比 并 更 新 集合 5S, 和 5; 。 每 隔 10 次 迭 


从 ， 


吧 


1 透镜 成 像 
Fig.1 Lens imaging 


集合 5 和 5; 分别 记 录 了 当前 已 知 的 适应 值 最 好 的 10 个 


解 和 其 反 向 解 , 而 集合 2 作为 % 和 3 记录 了 最 多 100 个 由 8 


和 5 的 解构 成 的 有 序 对 。 


较 小 的 "生成 的 镜像 范围 更 大 ， 而 较 大 的 9 
更 小 。 通 过 上 述 的 改进 之 


由 于 透镜 成 像 中 的 "是 控制 者 原生 解 和 镜像 之 间 的 距 


离 ， 
成 的 镜像 范围 


后 固定 的 n 值 已 经 很 难 适应 FPA 算 
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法 全 局 和 


局 部 搜索 随机 转换 的 机 制 ， 


趟 


体 公 


其 中 ，1 代表 当 


式 如 下 : 


常数 ， 控 制 
也 就 是 控制 


小 ， 


控制 常数 


得 测试 函数 ， 如 


线 变化 趋势 。 


的 取 


图 2 所 示 。 


通过 图 


就 被 放大 ， 


3 可 以 看 出 当 
n 作为 缩放 因 


子 碗 


原生 解 和 镜 
值 ， 选 取 Rastrigin’s 


人 民 初 期 ， 


n=(+(/T)2Y 
前 迭代 次 数 ，7 代表 最 大 迁 代 次 数 ， 为 控制 


李 大 海 


因此 需要 


改变 n 取 值 方 


(9) 


对 此 kk 值 能 
像 解 之 间 的 


3 次 对 比 ， 


够 控 币 


缩放 因子 的 大 


距离 ， 


所 以 为 了 精确 


函数 作为 对 比 不 同 值 


k 分 别 取 5、 
取 10 时 效果 最 佳 。 


10、15。 


n 值 被 放大 镜像 的 取 


算法 的 探索 


搜索 和 探索 速 


| ek 


的 镜像 取 值 范 转 
于 帮助 算法 跳出 
条 轮 和 上 一 轮 获 取 的 解 有 可 


也 变 小 ， 力 


局 部 最 优 解 。 


区 域 也 随 之 变 大 ， 
度 。 随 着 迭代 次 数 的 增 力 


0 强 算法 的 搜 


Nn 


有 利 
也 随 


索 精 度 ， 


人 


次 如 果 


会 一 


但 是 ; 


E 复 。 
一 样 的 ， 


也 就 


了 每 一 次 迭 


代为 一 下 


合 ,， 通过 把 每 次 迭代 9 


集合 中 ， 


个 体 。 


是 一 种 映射 


且 按 式 (8) 生 成 集合 
以 通过 


» 可 


通过 把 全 


mt 


能 
通过 自 适 应 收缩 因 
避免 了 重复 存储 
相 较 于 基本 的 透镜 成 像 策 
民生 成 当前 解 的 反 向 解 


寸 集合 2 中 的 二 
2 内 部 二 
与 反 向 解 作 对 比 ， 胜 出 者 作为 更 新 种 群 的 向 量 。 


的 次 端 。 
咯 ， 增 前 


0 .集合 


元 组 映射 的 


的 透镜 成 像 策 
透镜 成 像 策略 可 以 每 一 下 


各 的 每 


轮 生 成 
合生 成 100 个 有 


个 当前 


同 轮 的 反 向 解 找寻 到 最 合适 的 反 向 解 。 


策略 在 原 


有 的 其 


础 上 进 


多 样 性 。 


并 


增强 型 透镜 成 像 策 


更 好 的 适应 整个 搜索 过 程 。 


为 了 对 比 增强 


步 扩大 了 搜索 空 
且 对 比 透镜 成 像 策 略 的 
名 采取 自 


型 透镜 成 像 和 透镜 成 像 策略 


样 ， 罚 


模式 。 
E 成 的 解 和 反 向 解 存 入 5, 和 5， 


固定 缩放 攻 
适应 的 缩放 因子 ， 


算法 迭 


型 透镜 
而 是 


P 么 存储 在 8 
子 n 计算 不 同 轮 的 反 向 解 是 


于 算法 前 
之 变 小 ， 
同时 也 
代 到 后 期 ， 


成 像 策略 改变 
改 为 每 十 次 迭 


2 与 % 和 中 的 关系 就 


元 组 映射 到 5, 和 5; 内 的 


解 与 反 


向 解 和 当前 解 


解 的 反 
序 对 ， 
因此 增 
间 ， 


相 较 于 基本 
向 解 ， 增 强 型 
通过 搭配 不 
强 型 透镜 成 像 
有 利于 算法 的 


子 n， 本 文 的 
使 得 策略 可 以 


的 优 劣 性 ， 以 


基本 FPA 算法 为 基础 ， 选 取经 典 测试 函数 Sphere、Matyas、 
内 容 如 表 1 所 示 ， 适 


Rastrigin’s 和 Ackely’s 作为 参照 ， 详 细 


应 值 如 图 3 所 示 。 
表 1 测试 函数 
Tab. 1 Test function 
测试 函数 维度 ”范围 ”最 优 解 峰值 
f= 30  [-100,100] 0 ” 单 峰 
i=l 
f(x,y)=0.26*(x2 +y)—0.48* x*y 30  [-100,100] 0 ” 单 峰 
JOYy)=-20+ 妇 十 只 2 一 记 
10(cos2xx+cos2xy) 30 [5.12,5.12] 0 多 峰 
f (x,y)=—20exp(-0.2Y0.5(x2 + y7))— 30 [30,30] 0 多 蜂 
exp(0.5(cos(2xx) +cos(2zy)))+20+e ” 
Objective space 
10 T T 
Hk=10 
—e—k=5 
一 半 一 k=15 


Best score obtained so far 


1070 


10715 
0 


lteration 


图 2 对 


比 不 同上 值 结 


Fig.2 Compare the results of different k values 
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图 3 策略 适应 值 对 比 


Fig.3 Comparison of policy adaptation values 


从 图 3 中 可 以 很 明显 的 看 


和 收敛 速度 。 


2.2 基于 随机 邻 域 变异 策略 的 搜索 方式 

在 解决 复杂 问题 时 ， 
局 部 最 优 解 ，FPA 算法 往往 很 难 跳出 局 
方案 是 在 搜索 过 程 中 
041, 模拟 退火 [ 
E 和 


略 的 高 效 特 


上 ， 在 四 个 测试 函数 上 ， 增 强 
型 透镜 成 像 FPA 算法 比 透 镜 成 像 FPA 算法 有 更 好 的 收敛 性 


随 着 算法 进入 探索 后 期 ， 


旦 陷入 


尽量 增加 候选 解 


部 最 优 。 一 般 的 解决 
的 多 样 性 。 如 混沌 映射 


5 等 。 近 年， 因 差 分 进化 算法 


EDE)09 的 变异 策 


参数 简单 


已经 被 许多 学 者 | 


] 来 和 其 他 算法 结合 。 


DE 算法 


使 用 


“DE/best/1 ” 


这 意味 


DE/neighbor/1 


中 当 


中 最 好 的 候选 解 作为 3 
索性 相对 较 差 。Das08 进 
(DE/neighbor/1), 该 策略 引入 邻 域 的 信息 
敛 的 速度 和 收敛 的 精度 。 对 于 种 群 
体 中 选择 邻 域 集 
前 个 体 邻 域 中 的 最 好 一 个 作为 基 向 量 。 


[71。 


在 


有 更 好 的 探索 性 ， 但 开发 性 


的 基本 的 变异 策略 有 1 
前 者 策略 中 基 向 量 从 种 群 中 随机 选择 ， 


种 : “DE/rand/1 ”和 


较 差 。 后 者 则 获取 种 群 


基 向 量 ， 


步 


随机 的 从 群 


异 策略 可 以 表述 为 如 下 : 


革 中 9 Npest 


Vi= Xe + PF X(NK —X2) 


为 邻 域内 最 优 个 体 ; 


随机 挑选 的 个 体 。 


将 其 


其 中 ， gq e[0,1] 


应 用 于 FPA 算 的 


rand <q : x =xX*+rand(x; —xX!)+rand(xi 一 


优 解 ，xiw 是 邻 域内 最 优 解 。 


从 而 
是 出 了 


态 入 ,为 邻 域 内 剩余 个 体 中 


具有 较 好 的 开发 性 ， 但 探 
基于 邻 域 的 变异 策略 
指导 搜索 ,以 加 快 收 
中 的 每 个 个 体 在 每 一 代 ， 
合 ,， 并 以 基础 向 量 
DE/neighbor/1 变 


由 


(10) 


本 文 提出 了 一 种 基于 DE/neighbor/1 的 改进 变异 策略 , 并 


局 部 搜索 。 


rand > 9 = Xp + P(X —XI)+E (XY 


的 常数 ， 


x 对 代表 第 1+41 次 的 值 ， 
对 于 


改进 变异 策略 表述 如 下 : 
ax) GD) 

(12) 
x* 全 局 最 


式 (11) 中 汶 , 台 总 都 是 从 个 体 


一 发 ) 
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中 随机 选取 的 , 而 式 (12) 中 总 忆 , 必 都 是 从 邻 域内 随机 选取 的 个 


体 。 
对 于 邻 域 的 划分 以 当前 解 为 中 心 ， 


按 欧 氏 距离 由 小 到 大 


选取 邻 域 集合 ，n 维 的 欧 氏 距离 公式 如 下 所 示 。 
d= >》 (13) 
需要 注意 每 一 代 邻 域内 有 数量 限制 ， 而 算法 迭代 过 程 中 


的 二 选择 合理 的 邻 域 个 体 数量 N; ,对 平衡 全 局 和 局 部 搜索 能 
力 , 提高 算法 寻 优 能 力 起 着 重要 作用 。 受 文献 [19] 启 发 , 邻 域 
数量 按 如 下 公式 设置 : 


N; = Ni, + (Ni — Ni)x RAC )— fm m 
PF 7) fin) (14) 


中 : Ns 和 Ns 分别 代表 N 的 下 限 和 上 限 ， fi 代表 当前 种 群 


™ 


适应 值 的 最 小 值 ， f(%) 为 个 体 i 当 代 的 适应 值 。 

算法 初期 ， 种 群 中 优质 个 体 数量 较 少 ， 能 够 寻 得 全 局 最 
优 的 概率 较 小 ， 所 以 前 期 需要 加 强 全 局 搜索 能 力 ， 以 获得 较 
多 的 优质 个 体 , 提高 全 局 寻 优 的 能 力 。 根据 式 (14) 可 知 , 初期 


PFO) fn) 远大 于 f(%)-fis，N; 趋向 于 Ns 。 此 时 算法 具有 
较 好 的 全 局 搜索 能 力 ， 从 而 加 速 搜索 能 力 。 随 着 迭代 次 数 的 
增加 ， 种 群 中 优质 个 体 不 断 增加 ， 六 WCGo)- so 和 7CD- ja 


均 趋 向 于 0，N; 趋向 于 Ne， 此 时 该 算法 具有 较 好 的 局 部 开发 
能 力 ， 从 而 根据 邻 域 信息 寻找 最 优 解 。 

式 (10) 相 较 于 差分 进化 算法 中 的 单个 向 量 的 策略 ， 
两 个 差分 向 量 的 变异 策略 可 以 提高 种 群 的 多 样 性 ， 并 且 
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通过 自 适 应 邻 域 数 量 策略 ,动态 调整 种 群 邻 域 的 的 数量 
MSFPA 充分 利用 种 群 中 “随机 ”信息 和 邻 域内 “最 优 ” 信 息 ， 
合理 的 平衡 了 算法 的 全 局 和 局 部 搜索 能 力 ， 加 速 了 信息 的 交 
换 速 度 ， 增 加 种 群 的 多 样 性 ， 避 免 算法 陷入 局 部 最 优 解 的 同 
时 ， 加 速 了 算法 的 寻 优 速度 ， 而 且 没 有 引入 其 他 参数 ， 算 法 
的 复杂 度 没有 增加 。 


Sphere 函数 等 高 线 Matyas 函 数 等 高 线 
一 3 a pr 


-50 -一 一 - - -20 上 
50 40 30 20 -0 0 10 20 30 40 50 -20 -5 0 0 5 


图 4 策略 等 高 线 对 比 
Fig.4 Strategic contour comparison 
2.3 MSFPA 的 全 局 和 局 部 搜索 能 力 分 析 
算法 的 全 局 和 局 搜索 能 力 通常 作为 衡量 算法 性 能 的 两 个 
指标 。 全 局 搜索 代表 算法 的 前 中 期 的 收敛 能 力 ， 依 靠 全 局 搜 


能 够 
帮助 算法 跳出 局 部 最 优 解 。 而 且 使 用 全 局 最 优 解 作为 当代 的 
基 向 量 ， 有 利于 解 的 优良 基因 被 继承 。 而 式 (12) 中 的 代表 
变异 因子 ， 是 一 个 实 常数 因数 ， 控 制 差 分 向 量 缩放 的 作用 。 
对 于 算法 而 言 f 如 果 取 固定 值 的 话 太 小 或 者 太 大 都 会 导致 


索 算法 可 以 收敛 到 全 局 最 优 解 大 致 的 位 置 。 而 局 部 搜索 是 算 
法 后 期 提升 算法 的 搜索 精度 ， 能 够 让 算法 计算 出 的 解 无 限 接 
近 全 局 最 优 解 。FPA 算法 依靠 概率 转换 因子 P 去 切换 全 局 和 


算法 过 早 收敛 。 因 此 使 用 自 适应 变异 因子 Po 去 避免 不 合理 的 


参数 设置 带 来 的 影响 ， 有 具体 公式 如 下 : 
= fn f 
Fox = (Fins = Fin)X 人 
F= frr tin J (15) 
Fax trand x Pin i f 
其 中 :Fw 和 Fs 分别 为 F 的 上 下 限 ; rand 是 [0,1] 之 间 的 随机 


值 ，f; 为 个 体 坟 的 适应 度 值 ， 记 fi 为 当前 种 群 中 最 大 最 小 
的 适应 度 值 ， 了 为 当前 种 群 的 平均 适应 值 。 
式 (15) 在 文献 [20] 的 自 适 应 因子 的 基础 上 作出 改动 , 当 fi< 
7 时 不 在 是 固定 的 ,而 是 一 个 以 Fs 为 下 限 的 随机 值 ， 增 
加 算法 的 随机 性 。 自 适应 变异 因子 能 够 平衡 算法 寻 优 能 力 ， 
一 定 程度 上 降低 算法 对 参数 的 敏感 程度 ， 保 证 收敛 速度 和 避 
免 早熟 。 
为 了 进一步 验证 策略 的 有 效 性 ， 把 加 入 随机 邻 域 变异 策 


局 部 搜索 ,全 局 搜索 通过 Levy 飞行 增加 算法 的 随机 性 有 
利于 算法 的 前 期 前 中 期 探索 能 力 ， 最 后 通过 式 (1) 更 新 种 
群 , 有 利于 算法 朝 着 全 局 最 优 解 方向 探索 。 局 部 搜索 通过 
式 (3) 更 新 种 群 , 从 种 群 中 选择 最 好 的 候选 解 作为 基 向 量 ， 
保证 算法 的 局 部 搜索 阶段 具有 较 好 的 开发 性 。 而 MSFPA 
算法 引入 增强 型 透镜 成 像 ， 通 过 统计 10 轮 的 当前 解 以 及 
对 应 反 向 解 的 排列 组 合 后 选择 出 全 局 最 优 解 。 相 较 于 透 
镜 成 像 只 是 对 比 当 前 解 和 反 向 解 ， 增 强 型 透镜 成 像 策 略 
用 有 更 多 的 候选 解 供 其 选择 ， 增 加 算法 的 全 局 搜索 能 

以 及 种 群 的 多 样 性 , 并 且 通 过 对 缩放 因子 4 采取 自 适应 的 
方式 控制 当前 解 和 反 向 解 的 距离 ， 使 其 适应 整个 搜索 阶 
段 。 另外 增加 了 随机 令 域 变异 策略 , 在 花 授 粉 的 局 部 搜索 
阶段 ， 对 比 常 数 因子 94 和 rand 决定 何 种 方式 更 新 种 群 。 式 
(12) 因 为 使 用 全 局 最 优 解 作为 基 向 量 , 使 其 偏向 于 全 局 搜 


这 


略 的 FPA 和 标准 FPA 算法 进行 比较 ， 测 试 函数 如 表 1 所 示 ， 
对 比 图 如 4 所 示 。 


索 的 方式 ， 保 证 了 算法 即使 在 局 部 搜索 阶段 也 保持 一 
的 收敛 速度 ; 式 (13) 使 用 是 邻 域 集合 作为 基 向 量 , 并 且 集 


图 4 中 Sphere 和 Matyas 函数 分 别 选 取 了 第 100 代 和 第 
300 代 的 解 作为 对 比 ，Rastrigin’s 和 Ackely’s 选取 的 是 第 800 
和 第 1000 代 的 数据 作为 对 比 。 选 取 的 对 比 函 数 包含 单 峰 和 


合 内 部 的 数量 是 会 随 着 迭代 变化 更 好 的 适应 算法 前 期 收 
敛 速度 快 ， 后 期 收敛 精 度 高 的 需求 。 且 配 合 适应 的 变 
异 因 子 控制 算法 的 平衡 性 , 使 得 MSFPA 算法 有 着 更 好 的 


多 峰 ， 充 分 对 比 策略 对 于 不 同类 型 函数 的 适应 性 。 需 要 注意 
的 是 随机 邻 域 变异 FPA 在 Sphere 和 Matyas 函数 上 只 有 一 个 
点 ， 因 为 两 次 取 值 点 重合 了 ， 而 等 高 线 上 没有 FPA 的 点 是 
于 FPA 算法 获取 的 解 超出 了 范围 ， 也 就 是 说 随机 邻 域 变异 
FPA 在 第 100 次 迭代 的 时 候 已 经 收敛 到 0 附近 了 ， 而 同一 时 
期 的 FPA 算法 离 最 优点 很 远 。 对 于 Rastrigin’s 和 Ackely’s 玫 
数 而 言 分 别 取 了 第 800 代 和 第 1000 代 作 比较 ， 可 以 看 见 两 
次 取 值 都 是 不 如 随机 邻 域 变异 FPA 算法 。 


综合 寻 优 能 力 。 
2.4 MSFPA 算法 的 时 间 复 杂 度 分 析 

设 需要 优化 的 目标 函数 为 f(x) ， 解 空间 的 维度 为 4 ， 根 
据 FPA 算法 的 步骤 ，FPA 算法 的 时 间 复杂 s 度 为 OL+JCD) 。 
根据 MSFPA 算法 步 又, 种 群 规模 为 a, 初始 化 个 体位 置 的 时 
间 复 杂 度 为 0(rd) ， 根 据 初 始 位 置 生成 适应 值 的 时 间 为 f(4) 。 
全 局 搜索 中 基于 Levy 飞行 生成 步 长 的 时 间 为 引 ， 根 据 当 前 
位 置 产生 下 一 代 时 间 为 多; 局 部 搜索 中 的 变异 策略 产生 随机 
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的 时 间 为 6;。 


数量 级 另 


录用 定稿 


数 的 
变异 


时 间 为 所 ， 根 据 当 前 
因子 F 的 时 间 为 忆 ， 


| 
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科 置 产生 下 一 代 的 时 间 为 名 ,设置 
设置 邻 域 集合 的 时 间 为 点 。 


全 局 搜索 的 时 间 复 杂 度 为 
Olnd)+ On(f(d)+tad+%d+f(d))=0(d+f(d) (16) 

局 部 搜索 的 时 间 复 杂 度 为 
Ond)+On(f (D+Gtot+sts+f(d))=0d+f(d) (17) 


假设 得 到 新 的 适应 值 与 当前 适 


am， 经 过 比较 后 作为 当前 适应 值 的 时 
精英 集合 以 及 对 比 过 后 生成 新 适应 值 作 为 全 局 适应 值 


立 值 比较 的 时 间 为 
间 为 上 ， 


后 续 生成 


2*0O(d+f(d)) 
可 以 看 到 MSFPA 


| ， 增 强 型 透镜 成 像 以 及 随机 邻 
增加 算法 的 时 间 复 杂 度 。 


则 算法 总 的 时 间 复 杂 度 为 


+O(n(oi +0,+03)) =O(d 
与 FPA 算法 的 时 


改进 的 MSFPA 算法 伪 代 码 如 下 所 示 。 


算法 1 
输入 : 目标 函数 /Co ; 逢 
异 因子 下 
输出 : 全 局 最 优 解 。 


MSFPA 算法 
规模 NP ; 转换 概率 P 和 概率 因子 9 ; 变 
上 下 限 ( Fa, Fn )， 最 大 迭代 次 数 Maxlter 。 


+f(d) (18) 


间 复 杂 度 属于 同一 个 
或 的 搜索 策略 并 没有 


中 条 


初始 化 种 群 ， 在 初始 化 种 群 中 找到 全 局 最 优 解 x 


for t=1: Maxlter 
for i=1: NP 


if rand>P 
绘制 
按照 式 (1) 对 好 ” 


进行 全 局 授粉 阶段 
及 从 式 (2) 分 布 的 随机 步 长 向 量 


位 置 进行 更 新 


else 进行 局 部 授粉 


if rand<q 


在 当前 种 群 中 随机 选取 3 个 个 体 


按照 (12) 对 x 
else 
if 当前 解 的 适应 


行 更 新 


! 位 置 进 


值 >= 


mean( 适 应 值 ) 


按 式 (16) jf; >= 了 部 分 设置 变异 因 


else 


按 式 (16) 上 < 了 


end 


Ey 
2 


部 分 设置 变异 因子 忆 


按 式 (15) 设 置 邻 域 大 小 N 


o 


按 式 (14) 计 算 距 离 ， 选 取 离 当前 最 优 解 最 近 的 NN 个 解 
在 当前 种 群 中 随机 选取 3 个 个 体 


按 式 (13) 对 Xx! 位置 


end 


end 


行 更 新 


按 式 (18) 定 义 缩放 因子 


for i=1: NP 


按 式 (8) 生 成 当前 解 的 反 向 解 并 进行 越界 处 理 
对 比 当前 解 和 反 向 解 的 适应 值 ， 最 小 的 作为 


end 


if t<16 


把 当前 解 存 入 5 集合 
把 当前 解 的 反 向 解 存 入 5, 集合 


else if mod(t,16) 


~= 0 


当前 解 的 适应 值 


把 当前 解 和 当前 解 的 反 向 解 与 % 和 5, 中 的 解 对 比 ， 淘 汰 劣质 解 


else if mod(t,16) 


得 到 5S, 和 9 下 标 
for j = 1:size 


求 0 集合 二 元 组 对 应 % 和 5, 内 部 个 体 


的 有 序 对 集合 @ 
(2 ,1) 


者 适应 值 之 和 最 小 的 一 组 有 序 对 为 最 优 组 
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else 


第 一 元 素 映射 的 % 中 个 体 作 为 当前 解 


对 比 种 群 个 体 ， 


x = 最 优 个 体 ， 作 为 全 
3 ”实验 结果 和 分 析 


lll 
和 
改 
be 
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为 了 研究 多 策略 花 授 粉 算法 的 有 效 性 , 选取 CEC2013 测 


试 函 数 集 进 行 测试 ， 


依次 对 30/80 维度 进行 测试 ， 


各 算法 独 


立 运行 10 次 , 记录 最 优 值 、 平 均值 、 标 准 差 并 与 基本 的 花 授 
粉 算法 和 其 他 改进 算法 进行 对 比 。 参 与 对 比 的 算法 如 表 2 所 


示 ， 测试 函数 基本 信息 如 表 3 所 示 ， 适 应 度 | 


线 和 最 优 解 分 


布 如 图 5 和 6 所 示 。 测 试 函数 分 为 三 类 :Fl1-F5 单 模 态 函数 ， 


F5-F12 为 多 模 态 函数 ,F13-F15 为 复合 函数 。 


表 2 对 比 算法 

Tab.2 Comparison algorithm 

算法 算法 简称 
花 授粉 算法 FPA 

三 重 动态 调整 的 花 授粉 算法 HLFPA®?!] 
t- 分 布 扰动 策略 和 变异 策略 的 花 授 粉 算法 TMFPAD2] 
多 策略 融合 的 麻雀 算法 ISSADP3] 

表 3 测试 函数 CEC2013 说 明 


Tab. 3 Test function CEC2013 description 


Type Name 万 ff (x) 
Fl:Sphere Function -1400 
Unimodal Functions F2:Rotated Discus Function -1100 
F3:Different Powers Function -1000 
F4:Rotated Rosenbrock’s Function -900 
F5:Rotated Ackley’s Function -700 
F6:Rotated Weierstrass Function -600 
F7:Rastrigin’s Function -400 
F8:Roatetd Rastrigin’s Function -300 
Basic Multimodal . 
F9:Non-Continuous Rotated 
Functions a -200 
Rastrigin’s Function 
F10:Schwefel’sFunction -100 
Fll:Lunacek Bi_Rastrigin Function 300 
F12:Rotated Lunacek Bi_Rastrigin 2 
Function 
F13:Composition Function 2 
800 
(n=3,Unrotated) 
Composition F14:Composition Function 4 1000 
Functions (n=3,Unrotated) 
F15:Composition Function 8 
1400 
(n=3,Unrotated) 
注意 : 上 述 函 数 均 来 自 CEC2013 测试 集 

单 模 态 函数 Fl1-F5: 主要 特征 是 在 搜索 区 间 内 只 有 一 个 
局 部 最 小 值 。 

多 模 态 函 数 (F6-F12): 这 些 函 数 具有 多 个 局 部 最 优 值 ， 但 
是 与 每 个 函数 相对 应 的 单个 全 局 最 优 值 。 经 典 导数 的 优化 算 
法 在 此 类 函数 上 不 能 取得 满意 效果 ， 因 为 它们 会 陷入 局 部 最 
优 解 。 此 类 别 的 7 个 函数 有 以 下 特征 : 在 任何 地 方 都 是 不 可 


分 离 的、 不 对 称 、 连 续 的 ， 但 仅 在 一 组 点 上 是 可 区 


分 的 ; 在 


任何 地 方 都 是 连续 的 ， 但 在 任何 地 方 都 是 可 微分 的 


的 适应 值 并 记录 两 


if 最 优 组 第 一 元 素 的 适应 值 > 最 优 组 第 二 元 素 的 适应 值 


第 二 元 素 映 身 


的 $ 中 个 体 作 为 当前 解 


且 最 佳 和 次 佳 最 优 解 距离 较 远 。 
复合 函数 (F13-F15): 每 个 


合 函 数 均 


多 个 


局 部 最 优 到 全 局 最 优 的 狭窄 谷 ， 具 有 大 量 的 局 部 最 优 值 ， 并 


; 具有 从 


本 函数 组 


成 。3 个 复合 函数 都 是 多 模 态 ， 不 可 分 离 


的 和 部 队 称 的 ， 并 
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且 围绕 不 同 的 局 部 最 优 具 有 不 同 特征 。 
在 实验 中 ， 每 个 函数 都 有 一 个 相对 应 的 全 局 最 优 值 ， 该 
最 优 值 不 会 随 着 问题 的 维度 变化 。 实 验 所 有 函数 的 可 变 范围 
固定 为 [100,100]? ， 本 次 实验 对 15 个 函数 测试 每 种 算法 ， 并 
且 上 述 函 数 反复 随机 运行 10 次 。 

表 4、5 分 别 给 出 了 30 维 、80 维 的 实验 数据 。 从 表 4 中 
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融入 多 策略 后 该 算法 的 低 维度 问题 寻 优 能 力 得 到 了 提升 ， 并 
且 搜 索 精 度 方 面 得 到 了 加 强 。 

表 5 列 出 了 算法 在 80 维 测试 函数 上 的 测试 数据 。 基 
和 30 维 数据 相似 ， 除 了 F5 和 F10， 其 他 的 单 模 态 ， 多 模 态 
和 复合 函数 上 都 优 于 对 比 算法 ， 最 后 的 Ave Rank 值 为 1.33， 
排名 第 一 。 其 Ave Rank 值 和 30 维 接近 ， 说 明 维 度 提 升 并 没 


数据 可 以 看 出 F1-F15 的 函数 MSFPA 都 是 第 一 ， 这 意味 着 该 
算法 在 整个 搜索 过 程 中 寻 优 的 精度 和 稳定 性 都 要 好 于 对 比 函 
数 并 且 在 其 他 对 比 算法 搜索 已 经 停滞 的 时 候 MSFPA 算法 依 
旧 在 向 下 探索 说 明基 于 笛 卡 尔 积 的 透镜 成 像 策 略 和 随机 邻 域 
的 变异 策略 可 以 帮助 MSFPA 算法 跳出 局 部 最 优 解 。 在 F13- 
F15 复合 函数 上 均 是 排名 第 一 ， 这 说 明 MSFPA 算法 在 处 理 
复杂 问题 的 寻 优 能 力也 强 于 对 比 函 数 并 且 对 比 其 他 算法 
MSFPA 算法 的 收敛 速度 和 收敛 精度 都 要 更 高 。 在 30 维 测 试 
情况 下 ，MSFPA 算法 的 Ave Rank 为 1， 排 名 第 一 。 这 表明 


一 一 pa Liat 
Mp i 
“hllpa Mi 
一 一 SSA 6 一 一 wp 
一 MSFPA i 
a 一 wsFPA 


1 


人 AAAAAAAAAAGAA 


Fl F2 


有 降低 该 算法 的 性 能 。 

5 是 在 30 维 的 情况 下 , 各 类 算法 的 收敛 情况 , 可 以 明 
显 看 出 ， 列 出 的 图 集中 在 单 模 态 ， 多 模 态 和 复合 函数 上 ， 殿 
收敛 精度 都 要 优 于 大 部 分 的 改进 算法 .MSFPA 利用 透镜 成 像 
和 随机 邻 域 的 变异 策略 扩展 搜索 范围 和 改进 搜索 精度 。 例 如 
F5、F6、F7、F8、F9 和 F10， 其 他 函数 都 基本 停滞 的 时 候 ， 
但 是 MSFPA 算法 还 在 向 下 搜索 。 这 表明 改进 后 确实 对 算法 
的 寻 优 精度 有 显著 提升 ， 而 且 该 算法 的 收敛 速度 对 比 其 他 的 
算法 也 相对 较 快 。 


一 一 和 IE 本 
要 Mip: 本 Mfpa 
一 一 wp 25| 一 一 ne 
和 ISSA A ISSA 
MsFPA | 一 wsFp 
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图 5 适应 度 曲线 


Fig.5 Adaptability curves 


图 6 是 所 有 算法 1000 次 迭代 获取 最 优 解 的 箱 线 图 ， 
Fl、F6、F8、F10 和 F13 中 的 MSFPA 箱 体 较 宽 ， 代 表 该 
算法 的 所 有 最 优 值 波动 情况 大 ， 也 就 是 算法 收敛 速度 较 
快 导 致 每 一 代 的 最 优 解 之 间 跨 度 较 大 。 其 他 的 测试 函数 
中 箱 体 比较 窗 ， 代 表 算 法 从 开始 搜索 到 迭代 结束 获取 的 


所 有 解 变 化 不 大 。 从 图 7 可 以 看 出 MSFPA 在 F1、F2、 
F6、F7、F8、F1l1、F13、F15 中 箱 体 的 下 限 比 其 他 算法 更 
低 ， 也 就 代表 它 的 搜索 精度 更 高 。 
在 15 个 测试 函数 中 ，30 维 有 15 个 ，80 维 有 13 个 测试 
函数 是 优 于 FPA、HLFPA、TMFPA、ISSA。 这 可 以 说 明 增强 
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3.1 Friedman 


为 了 验证 MSFPA 算 济 
24 对 算法 
种 非 参数 双向 方差 检验 方法 ， 


Friedman 检验 


检验 


李 大 海 ， 


型 透镜 成 像 策略 和 随机 邻 域 的 变异 策略 有 助 于 跳出 
解 ， 并 能 够 指导 算法 后 续 的 探索 。 


局 部 最 优 


和 其 


他 算法 的 显著 性 差异 ， 采 用 


FE 进行 参数 检验 。Friedman 检验 是 一 
一 般 用 来 检测 数据 之 间 的 显 


将 检验 MSFPA 算法 和 对 比 算 法 在 30 维 
检验 结果 如 表 6 所 示 。 表 


等 : 多 策略 增强 花 授 粉 算法 及 其 应 用 


法 之 间 的 差异 性 就 越 大 ， 当 
是 服从 k-1 自由 度 的 分 布 。 
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n 和 k 足够 大 时 (一 般 n>10,k>5)， 


127 2 5 
F, RD k(k—1)?/4] 


(19) 


，80 维 


的 ; 


Pp 的 P-value 表示 渐进 显著 性 ， 


著 性 差异 ， 多 用 来 做 单 目标 之 间 的 优 劣 比较 。 一 般 实现 步 。 ”判断 是 否 存在 显著 性 差异 的 重要 指标 , 如 果 小 于 0.01, 则 表示 
又 如 下 : 各 项 数据 之 间 存 在 显著 性 差异 。 表 中 其 他 值 表示 各 算法 的 秩 平 
a) 收集 实验 数据 。 列 出 每 个 问题 i 的 最 好 与 最 差 结果 排 ” 均值 。 表 中 的 30、80 维 的 P-value 值 明 显 小 于 0.01。 这 表明 
名 ，mQdxK< 方 。 MSFPA 算法 和 其 他 对 比 算 法 存在 显著 性 差异 ， 这 种 差异 有 可 
b) 求 出 所 有 算法 的 平均 排名 ， 并 计算 出 最 终 排名 。 能 是 多 策略 带 来 的 解 的 多 样 性 以 及 搜索 能 力 的 加 强 。 各 算法 的 
&=d-0Z257 。 秩 平均 值 越 小 则 代表 算法 性 能 越 好 。 秩 平均 值 中 , MSFPA 在 每 个 维度 都 是 最 小 的 , 因此 MSFPA 算 
Friedman 统计 值 计算 公式 如 (12) 所 示 ，# 的 值 越 小 , 各 算 。 ”法 在 各 维度 上 的 优化 能 力 都 得 到 了 提升 。 
， : 。 "| : 
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图 6 最 优 解 箱 线 图 
Fig.6 Optimal solution boxplot 
4 ”汽车 动力 传动 参数 优化 a) 动力 性 能 发 挥 程度 的 评价 指标 为 驱动 功率 损失 率 。 
人 到 b) 经 济 性 能 发 挥 程度 的 评价 指标 为 有 效 效率 利用 率 。 
由 内 燃 机 理论 可 知 ,汽车 的 动力 性 和 燃油 经 济 性 指标 是 相 c) 汽车 动力 传动 系统 匹配 的 综合 指标 为 汽车 能 量 利用 率 。 
互 矛盾 的 。 即 无 法 脱离 汽车 的 动力 性 而 单纯 追求 燃油 经 济 性 ， 优化 模型 的 设计 变量 选 为 
也 无 法 摆脱 燃油 经 济 性 追求 动力 性 , 因此 只 能 在 汽车 的 动力 性 区 =[X Xs, Xs, Xa Xs Xo = Liisa io ionics i (20) 
和 燃油 经 济 性 选择 一 个 平衡 方案 .而 汽车 动力 性 燃油 经 济 性 的 。 ”其 中 : is 为 变速 器 第 j 档 的 传动 比 (j=42.…5); in e[5] ,za el 入， 
综合 评价 体系 和 指标 包括 : ia E[0.5,3] ,is e[0.5,2] , is s[034] ,加 e[2,6] 为 主 减速 器 传动 比 。 
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表 4 实验 结果 


Tab. 4 Experimental results 


测试 函数 Mean/Std/Rank FPA TMFPA HLFPA ISSA MSFPA 
fl Mean/Std 1.60E+03/2.68E+02 1.00E-01/9.50E-02 8.58E+03/1.78E+03 9.93E+00/1.88E+01 3.30E-12/1.91E-12 
Rank 4 2 5 3 1 
f2 Mean/Std 9.24E+04/2.07E+04 3.13E+03/1.66E+03 7.78E+04/1.65E+04 3.76E+04/5.23E+03 1.32E+02/1.25E+02 
Rank 3 用 4 有 1 
Mean/Std 1.42E+02/4.18E+01 7.15E-01/3.46E-01 2.24E+03/3.25E+02 1.80E+02/5.37E+01 5.93E-08/4.15E-08 
f3 
Rank 3 2 5 4 1 
f4 Mean/Std 1.36E+02/2.79E+01 4.34E+01/2.77E+01 5.71E+02/1.27E+02 1.48E+02/3.56E+01 2.01E+01/1.88E+01 
Rank 3 2 5 4 1 
f5 Mean/Std 2.10E+01/6.17E-02 2.10E+01/6.63E-02 2.10E+01/4.96E-02 2.10E+01/4.45E-02 2.10E+01/7.18E-02 
Rank 3 用 4 5 1 
f6 Mean/Std 3.90E+01/1.53E+00 2.58E+01/5.09E+00 3.95E+01/1.99E+00 4.40E+01/2.71E+00 2.21E+01/4.25E+00 
Rank 3 2 4 本 1 
f7 Mean/Std 2.55E+02/2.84E+01 1.09E+02/2.70E+01 3.54E+02/2.95E+0 8.40E+02/7.62E+0 4.29E+01/1.39E+01 
Rank 3 大 4 3 1 
f8 Mean/Std 2.96E+02/2.21E+01 1.97E+02/3.95E+01 3.87E+02/2.92E+0 9.38E+02/2.41E+0 6.94E+01/2.94E+01 
Rank 3 也 4 5 1 
Mean/Std 3.16E+02/2.54E+01 2.44E+02/3.47E+01 3.71E+02/2.02E+0 1.03E+03/8.36E+0 1.27E+02/2.75E+01 
f9 
Rank 3 2 4 入 1 
Mean/Std 4.38E+03/6.58E+02 4.14E+03/5.15E+02 7.23E+03/2.82E+02 5.86E+03/7.56E+02 2.89E+03/5.20E+02 
f10 
Rank 3 2 5 4 1 
f11 Mean/Std 4.01E+02/5.63E+01 2.47E+02/5.85E+01 5.70E+02/3.49E+0 8.35E+02/1.89E+0 1.23E+02/1.52E+01 
Rank 3 2 4 5 1 
f12 Mean/Std 3.94E+02/4.29E+01 3.03E+02/4.46E+01 5.77E+02/5.26E+0 8.16E+02/3.97E+0 1.47E+02/3.86E+01 
Rank 3 用 4 3 1 
f13 Mean/Std 5.20E+03/9.38E+02 4.56E+03/3.92E+02 7.49E+03/4.03E+02 7.64E+03/7.81E+02 3.73E+03/7.74E+02 
Rank 3 2 4 5 1 
f14 Mean/Std 2.81E+02/6.23E+00 2.72E+02/7.28E+00 2.95E+02/5.77E+00 7.53E+02/4.09E+0 2.64E+02/7.69E+00 
Rank 3 2 4 有 1 
Mean/Std 1.79E+03/1.42E+02 1.09E+03/1.02E+03 2.35E+03/1.86E+02 7.71E+03/1.08E+03 4.19E+02/3.75E+02 
f15 
Rank 3 2 4 六 1 
Count 0 0 0 0 15 
Ave Rank 32 受 4.2 4.53 1 
Total Rank 3 2 4 5 1 
维度 :30， 以 上 数据 独立 运行 10 次 的 均值 
表 5 实验 结 
Tab.5 Experimental results 
测试 函数 Mean/Std/Rank FPA TMFPA HLFPA ISSA MSFPA 
fl Mean/Std 5.70E+04/1.13E+04 1.18E+03/4.25E+02 9.70E+04/1.15E+04 1.97E+04/6.95E+03 S.10E-02/1.68E-02 
Rank 4 2 5 3 1 
Mean/Std 1.34E+05/2.47E+04 3.63E+03/1.47E+03 3.08E+05/2.79E+04 1.24E+05/1.10E+04 9.99E+02/9.58E+02 
f2 
Rank 4 2 5 3 1 
f3 Mean/Std 8.06E+03/1.30E+03 5.02E+02/3.93E+02 1.98E+04/3.36E+03 2.40E+03/5.37E+02 4.26E-01/1.38E-01 
Rank 4 2 5 3 1 
f4 Mean/Std 6.05E+03/1.50E+03 5.13E+02/8.41E+01 1.37E+04/3.43E+03 2.26E+03/5.99E+02 2.53E+02/5.81E+01 
Rank 4 2 5 3 1 
/5 Mean/Std 2.13E+01/2.49E-02 2.13E+01/2.96E-02 2.13E+01/3.04E-02 2.11E+01/3.64E-02 2.13E+01/3.60E-02 
Rank 4 3 1 5 
Mean/Std 1.09E+02/8.19E+00 9.23E+01/8.74E+00 1.20E+02/5.52E+00 1.32E+02/2.84E+00 5.71E+01/1.13E+01 
f6 
Rank 3 2 4 5 1 
17 Mean/Std 1.64E+03/1.28E+02 8.22E+02/1.05E+02 1.91E+03/1.05E+02 1.96E+03/6.65E+01 3.56E+02/3.42E+01 
Rank 3 2 4 5 1 
Mean/Std 1.62E+03/1.84E+02 1.04E+03/1.24E+02 1.81E+03/1.10E+02 2.05E+03/3.20E+01 4.22E+02/8.46E+01 
f3 
Rank 3 2 4 5 1 
Mean/Std 1.69E+03/1.60E+02 1.16E+03/1.44E+02 1.82E+03/1.84E+02 2.39E+03/7.20E+01 5.63E+02/9.94E+01 
f9 
Rank 3 2 4 5 1 
A10 Mean/Std 2.08E+04/6.63E+02 1.72E+04/1.66E+03 2.43E+04/5.29E+02 1.47E+04/2.72E+03 1.69E+04/9.31E+02 
Rank 4 3 5 1 2 
Mean/Std 4.21E+03/7.29E+02 1.44E+03/1.77E+02 2.53E+03/8.46E+01 2.25E+03/3.23E+01 7.3SE+02/6.6SE+01 
f1l 
Rank 5 2 4 3 1 
Mean/Std 3.69E+03/5.46E+02 1.71E+03/2.24E+02 2.51E+03/8.96E+01 2.24E+03/2.85E+01 7.85E+02/5.49E+01 
f12 
Rank 5 2 4 3 1 
f13 Mean/Std 2.15E+04/7.56E+02 1.83E+04/1.04E+03 2.49E+04/5.49E+02 2.06E+04/1.16E+03 1.71E+04/1.47E+03 
Rank 4 2 5 3 1 
Mean/Std 4.46E+02/1.87E+01 4.86E+02/1.77E+01 4.69E+02/9.28E+00 3.00E+03/5.00E+02 3.79E+02/2.19E+01 
f14 
Rank 2 4 3 5 1 
Mean/Std 7.64E+03/4.64E+02 6.72E+03/2.04E+03 9.09E+03/6.57E+02 2.12E+04/3.61E+02 2.21E+03/2.06E+02 
f15 
Rank 3 2 4 5 1 
Count 0 0 0 13 
Ave Rank 3.67 2.2 4.27 3.53 1.33 
Total Rank 4 2 4 肝 1 


维度 :80， 以 上 数据 独立 运行 10 次 的 均值 
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4.1 目标 函数 其 中 ，Tons =132， 发 动机 最 大 转 矩 IN*m); fzo = ， 驱 动 轮 上 
目标 函数 为 的 法 向 反作用 力 (N); 2 =0.7 为 地 面 附着 系数 。 
F(X)=Af 0 +h CD) ge(X)=0.85xg—in /i <0 (29) 
其 中 : 4 为 动力 性 发 挥 程度 加 权 因 子 ; 4 为 经 济 性 加 权 因 子 ; 其 中 ， 公 比 y=s 语 人 。 其 他 约束 如 下 所 示 。 
上 CD) 为 动力 性 分 目标 函数 ， 记 (» 为 经 济 性 分 目标 函数 。 选 择 87(X)=isi /iss —1.20g <0 (30) 
0<4<%hb<1， 取 加 权 因 子 41414=0.2， 罗 =0.8 。 gs(X)=0.80g is /isa <0 (G31) 
动力 性 分 目标 函数 : go(X)=is /iss -1.10g <0 (32) 
gio(X)=0.75g is /i <0 (33) 
f= d= CE _F 3 (22) gu(X)=ia /is -1.05g<0 (34) 
56=1.06+0.04i2? ， 五 万, 及 分 别 为 汽车 的 驱动 力 、 滚 动 阻力 和 空 gn(X)=0.70g -isa /iss £0 (35) 
气 阻 为 5 ga(X)=iss/is -1.0g<0 (36) 
经 济 性 分 目标 函数 (在 速度 wu 下 行驶 某 段 距离 AS 的 耗 油 量 ): gu(X)=isz iss -0.9 /is 0 (37) 
KxPxg.x(n,P), gis(X)=is3/iss —0.95is /iss <0 (38) 
名 102xu, xD 和 (23) gie(X)=ica /iss —0.95is /i <0 (39) 
其 中 ，P 为 燃油 重度 ( N/L)， 取 值 为 7.0。K 为 加 权 系 数 ， 等 gn(X)=X,—-X<0 (40) 
速 时 取 1， 加 速 时 取 1.05; wu e[10,20,30,40,50] 为 汽车 行驶 加 速 gis(X)=X-X, <0 (41) 
(km/h)。 发 动机 的 燃油 消耗 率 变 化 如 图 7 所 示 。 其 中 ，X, 和 X, 为 变速 器 和 主 减 速 器 传动 比 上 下 限 值 构成 的 向 量 。 
表 6 Friedman 检验 结 4.3 汽车 传动 参数 优化 实验 结果 分 析 
Tab.6 Friedman test results 在 算法 性 能 的 基础 上 ， 对 多 策略 的 花 授 粉 算法 (MSFPA) 
维度 (D) P-value FPA TMFPA HLFPA ISSA MSFPA 在 汽车 动力 传动 参数 优化 问题 ， 结 合 之 前 的 测试 算法 进行 综 
30D 5.3308E-11 3.2 2 42 4.53 1 合 对 比 。 以 此 表明 该 算法 在 实际 应 用 问题 上 的 可 适用 性 。 具 
80D 5.8656E-7 3.67 2.2 427 3.53 1.33 体 的 目标 函数 和 约束 条 件 已 在 第 二 节 中 有 详细 描述 。 对 比 结 
等 燃油 消耗 曲线 和 等 功率 曲线 的 三 维 拟 合 图 0 万 有 特性 曲线 果 如 图 8 所 示 ， 具 体 优 化 后 的 值 见 表 7。 
500 127 ek 
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200 地 
§ 100 E ® 
CC ， | 
600 
四 WU 
400 
2000 
200 00 = 0 
Ttq(N*m) 1000 n(r/min) 1500 2000 2500 20 40 60 80 100 120 140 
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图 7 燃油 消耗 率 变化 图 8 汽车 传动 参数 各 函数 适应 值 对 比 


Fig.7 Change of fuel consumption rate 


4.2 参数 优化 约束 条 件 


transmission parameters 


Fig.8 Comparison of adaptive values of functions of automobile 


在 燃油 汽车 优化 满足 动力 性 需要 , 之 后 满足 燃油 经 济 性 9: 图 8 为 了 方便 展示 放大 了 前 70 次 迭代 部 分 ， 从 图 中 可 
1X)=U 一 om SO (24) ”以 看 出 ，MSFPA 算法 的 收敛 速度 最 快 , 仅 在 第 5 次 迭代 左右 
其 中 ，w e[0,50] ， 汽 车 行驶 车 速 (km/b)。 就 已 经 到 达 最 低 点 ， 而 fpa 算法 虽然 初始 时 是 最 低 的 ， 但 是 
ga(X)=i i <0 (25) ”后 续 的 收敛 速度 明显 很 慢 。 而 其 他 的 算法 不 管 是 收敛 精度 还 
其 中 ， ee 是 在 收敛 速度 上 都 逊色 于 MSFPA 算法 。 
g3(X)=D, -Dw <0 (20) 根据 表 7 优化 后 参数 的 结果 可 以 看 出 ，5 种 对 比 算法 最 
其 中 ， DD 为 汽车 最 大 动力 因 es 后 适应 值 最 小 的 是 MSFPA 为 374， 其 他 算法 略 大 于 该 算法 。 
ga(X)=D 一 Do <0 (27) 而 汽车 传动 系 的 变速 箱 的 速 比 和 主 减速 器 的 速 比 X~Xs。， 其 
其 中 ， D, 为 汽车 最 高 挡 的 最 大 动力 因数 要 求 的 下 限 值 。 值 都 优 于 对 比 算法 。 因 此 MSFPA 算法 在 汽车 动力 传动 系统 
ee (28) 参数 优化 方面 相 较 于 对 比 算法 有 一 定 的 优越 性 ， 在 这 一 半 问 
7 题 上 具有 参考 价值 。 
表 7 ”优化 参数 详情 
Tab.7 Optimization parameters details 
算法 X, X, 5 Xs Xs Xs Fn 
FPA 1.28E+00 1.00E+00 8.34E-01 7.44E-01 7.31E-01 2.16E+00 3.76E+02 
MSFPA 4.54E-01 3.70E-01 3.22E-01 3.02E-01 3.02E-01 7.27E-01 3.74E+02 
TMFPA 1.56E+00 1.00E+00 6.78E-01 5.00E-01 4.43E-01 2.11E+00 3.75E+02 
HLFPA 1.28E+00 1.03E+00 9.23E-01 8.77E-01 8.77E-01 2.35E+00 3.77E+02 
ISSA 1.53E+00 1.00E+00 6.89E-01 5.00E-01 4.56E-01 3.00E+00 3.75E+02 
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为 了 改善 MSFPA 算法 的 收敛 速度 ， 收 敛 精度 以 及 容易 
陷入 局 部 最 优 解 的 问题 ， 提 出 一 种 结合 笛 卡 尔 积 透镜 成 像 和 
随机 邻 域 变异 策略 的 花 授粉 算法 ， 利 用 透镜 成 像 策 略 增加 解 
的 多 样 性 帮助 算法 跳出 局 部 最 优 解 ， 引 入 随机 邻 域 变异 策略 
加 强 算 法 的 局 部 搜索 能 力 ， 最 后 根据 实验 结果 对 比 4 种 不 同 
的 算法 ，MSFPA 算法 对 各 类 单 模 态 、 多 模 态 、 复 合 、 低 维 、 
高 维 都 具有 更 好 的 搜索 精度 和 收敛 速度 ， 最 后 结合 汽车 传动 
参数 优化 问题 进一步 验证 MSFPA 算法 的 实际 应 用 能 力 ， 结 
果 表 明 该 算法 具有 一 定 的 研究 意义 ， 为 求解 复杂 问题 提供 了 
种 新 的 思路 和 解决 办 法 。 
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